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Abstract 

Global and regional land carbon storage has been significantly affected by increasing atmospheric CO2 concentration and 
clinnate change. Based on fully coupled climate-carbon-cycle simulations from the Coupled Model Intercomparison Project 
Phase 5 (CMIPS), we investigate sensitivities of land carbon storage to rising atmospheric CO2 concentration and climate 
change over the world and 21 regions during the 130 years. Overall, the simulations suggest that consistently spatial 
positive effects of the increasing CO2 concentrations on land carbon storage are expressed with a multi-model averaged 
value of 1.04PgC per ppm. The stronger positive values are mainly located in the broad areas of temperate and tropical 
forest, especially in Amazon basin and western Africa. However, large heterogeneity distributed for sensitivities of land 
carbon storage to climate change. Climate change causes decrease in land carbon storage in most tropics and the Southern 
Hemisphere. In these regions, decrease in soil moisture (MRSO) and enhanced drought somewhat contribute to such a 
decrease accompanied with rising temperature. Conversely, an increase in land carbon storage has been observed in high 
latitude and altitude regions (e.g., northern Asia and Tibet). The model simulations also suggest that global negative 
impacts of climate change on land carbon storage are predominantly attributed to decrease in land carbon storage in 
tropics. Although current warming can lead to an increase in land storage of high latitudes of Northern Hemisphere due to 
elevated vegetation growth, a risk of exacerbated future climate change may be induced due to release of carbon from 
tropics. 
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Introduction 

Variations in the whole terrestrial carbon cycle are accompa- 
nied with increasing atmospheric CO2 concentration and climate 
change. The physical climate can influence the terrestrial carbon 
storage and the carbon exchange between atmosphere and land 
[1-5], and the sequent changes in atmospheric concentration of 
CO2 simultaneously affect the physical climate system [1,6]. By 
the end of the twenty-first century, there is an additional CO2 
change between 20 and 200 ppm considering warming alone, 
while the higher CO2 values can lead to an additional climate 
warming ranging between 0.1° and 1.5°C [1]. Zeng et al. (2004) 
[7] also suggested a positive feedback to global warming from the 
interactive carbon cycle introduces an additional increase of 
90 ppm in the atmospheric CO2 and then 0.6 degree additional 
warming during the period of 1860-2100 with the prescribed 
IPCC-SRES-AIB emission scenario. Regionally, using the Cou- 
pled Climate Carbon Cycle Model Intercomparison Project 
(C MIP), Cox et al. (2013) [8] estimates that warming alone will 
release 53±17 gigatonnes K~ of carbon to the atmosphere over 



tropical land from 30° north to 30° south. On the other side, 
increasing CO2 concentrations would influence land carbon 
balance, such as the stimulation of carbon storage, increases in 
gross primary production, net primary productivity and hetero- 
trophic respiration [6,8—12]. With increasing CO2 concentrations 
alone there is a widely distributed terrestrial carbon sink of 1.4— 
3.8PgCyr"' during the 1990s, rising to 3.7-8.6 PgCyr"' a 
century later [13]. Such an enhancement of land carbon storage 
was accessed due to CO2 fertilization effects with a rate of 
0.07PgCyr~ using simulation from Organizing Carbon and 
Hydrology considering fire disturbance [12]. 

As sensitivities of land carbon storage are largely disturbed and 
altered by both increasing CO2 concentrations and climate change 
[12,14—18], previous studies have focused on it [1,14]. Generally, 
when climate change is only accounted for, a reduction is 
suggested for the efficiency of the terrestrial ecosystems to absorb 
anthropogenic carbon emissions [2,19,20] and general negative 
sensitivities of land carbon storage have been revealed from the 
earth system models (ESMs)[14] and Dynamic Global Vegetation 
Models (DGVMs) [19]. Arora et al. (2013) [14] estimated earth 
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system models (ESMs) simulations from CMIP5, which represent 
the interactions between the carbon cycle and physical climate 
system. Negative impacts of continuous warming on terrestrial 
carbon storage have been suggested. Sitch et al. (2008) [19] used 
the five Dynamic Global Vegetation Models (DGVMs) and IPSL 
GCM to evaluate the terrestrial carbon cycle and climate-carbon 
cycle feedbacks, and found the similar sign of sensitivities of 
terrestrial land carbon to the changes in atmospheric CO2 and 
climate. Generally, terrestrial carbon balance and storage is 
discovered to be sensitive to rising atmospheric CO2 concentration 
and climate change [17,21,22]. Land carbon sensitivities to rising 
CO2 are positive implying increases under enhanced atmospheric 
CO2 concentration about 1.19~1.32 PgC ppm~\ while its 
sensitivities to temperature are negative ranging between 

-137 85 PgCK">[23]. 

Although 1)()th the common results suggest the positive response 
of terrestrial carbon storage to increasing CO2 concentration and 
its negative response to climate change [1,8,23], many potentially 
uncertainties are remained [14,19], such as the response 
magnitudes and regional pattern, often because they are poorly 
understood [24], or scantily quantified at scales relevant for 
models [1] . Several previous studies have been carried in order to 
evaluate and predict land carbon sensitivities to CO2 and 
temperature [1,14,23]. But most of them have not distinguished 
the regional variability due to impacts of rising atmospheric CO2 
concentration or climate change on carbon storage in order to 
identify the sensitive regions, respectively. On the basis of this 
point, we focus on the sensitivities of land carbon storage to the 
rising atmospheric CO2 concentrations and climate change at 
global and regional scales by multi-models, and compare the 
simulated variations between models that represent a potential 
source of uncertainty. Among regions which are more fragile to 
climate change, we quantify sensitivities of land carbon storage at 
regional scales. It is important for mitigation and adaptation of 
terrestrial ecosystems to climate change. In this study, we use 
simulations from earth system models offered by the CMIP5 to 
quantify and compare the carbon sensitivities to changes in CO2 
and climate for the period from 1860 to 1989. The disturbances of 
none-C02 greenhouse gases, land use and fire are totally ignored. 

Experiments and Methods 

The fifth phase of Coupled Model Intercomparison Project 
(CMIP5) can provide a common framework [25] to compare and 
assess land carbon storage responses in the context of climate 
simulations [14]. There are two experiments used in this analysis, 
which have been downloaded from http://cmiplOOpcmdi.llnl. 
gov/cmip5/forcing.html. These chosen experiments from CMIP5 
have been described by Peng et al.(2013; 2014) [26,27]. They are 
detailed as follows: (1) in the experiment of only considering the 
single effect of atmospheric increasing CO2 concentrations, 
biogeochemistry only responds to tlu- increasing CO2 concentra- 
tions in land models while the radiative forcing is fixed at the pre- 
industrial values in the atmospheric modules; (2) in the experiment 
of only thinking of the single effect of climate change, the radiative 
forcing responds to the increasing atmospheric CO2 concentration 
(prescribed by 1% increase of CO2 concentration per year in 
atmosphere from the pre-industrial values to the quadruple) while 
the biogeochemistry remains at the pre-industrial values in the 
biogeochemistry modules. In addition, in order to estimate and 
compare the sensitivities of land carbon storage to increasing CO2 
concentration and climate change, we use the simulations from six 
of the fully carbon-climate coupled ESMs, which participate in the 
CMIP5 intercomparison project. These models include Had- 



GEM2-ES [28], IPSL-CM5A-LR [29], CESMl-BGC [14], MPI- 
ESM-LR [30,31], CanESM2 [32] and BCC-CSMl-1 [33]. 

The ESMs, which are used to analyze the sensitivities, except 
for HadGEM2-ES and MPI-ESM-LR, do not include dynamical 
vegetation cover. Their spatial distribution is controlled by the 
competition het\\een different PFTs (Plant functional types) [14]. 
A patch-based representation of vegetation structure competition 
has been made [34]. HadGEM2-ES can simulate the transient 
shifts and geographic patterns of vegetation. Its atmosphere/land 
component resolution is 1.6° by 1.6°. A mixture of 12 PFTs is 
shown in HadGEM2-ES [28]. Over each patch a transient 
dynamics of vegetation also has be presented by MPI-ESM-LR. It 
represents the resolution of atmospheric component of 1.9° by 1.9° 
[31]. 

The IPSL-CM5A-LR [29] is the new generation Earth System 
Model developed by the Institute Pierre Simon Laplace. A 
resolution of 3.6° by 1.8° at latitudes and longitudes is used in the 
land and atmospheric components. There are 39 vertical levels of 
atmosphere. With a daily time step land component ORCHIDEE 
[35] simulates processes of photosynthesis, carbon allocation, litter 
decomposition, soil carbon dynamics, maintenance and growth 
respiration and phenology for 13 different plant functional types 
[34]. 

Version 1 of the Community Earth System Model (CESMl) is 
the successor to version 4 of the Community Climate System 
Model (CCSM4). In this model, global climate model consisting of 
land, atmosphere, ocean, and sea-ice components has been fuUy 
coupled [36]. The terrestrial carbon cycle is represented by the 
CLM4 land surface scheme. And it contains 16 different PFTs. 
Coupled carbon-nitrogen dynamics are included and expressed as 
CLM4CN [37,38]. 

CanESM2 is the new generation Earth System Model 
developed from the first generation Canadian earth system model 
(CanESMl) [21,39]. It is produced and developed by the 
Canadian Centre for Climate Modelling and Analysis (CCCma) 
and described by Arora et al. (2013) [14]. The atmospheric 
component has horizontal resolution of about 2.8° in a grid cell. 
The Canadian Terrestrial Ecosystem Model (CTEM) [40] is used 
to model terrestrial ecosystem processes. It simulates terrestrial 
carbon supported by three live vegetation pools (leaves, stem, and 
root) and two dead pools (litter and soil organic carbon) for PFTs 
(e.g. needle leaf evergreen and deciduous trees, broadleaf 
evergreen and cold and dry deciduous trees, and C3 and C4 
crops and grasses) [14]. 

BCC-CSMl.l is a fully coupled global climate-carbon model 
including interactive vegetation and global carbon cycle [33] . The 
atmospheric comporu-nt BCC-AGCM2.1 is a global spectral 
model with a horizontal resolution of 2.81 degrees, with vertical 26 
levels of atmosphere. For the land surface processes, the land 
surface component Atmosphere-Vegetation Interaction Model 
(A VIM) [41,42] is incorporated into the biogeophysical frame of 
NCAR/CLM3. With 15 PFTs including natural vegetation and 
crop, each a grid cell can contains up to four PFTs. Terrestrial 
carbon cycle components simulates biochemical and physiological 
processes. For example, photosynthesis and respiration of vegeta- 
tion, allocation of carbohydrate to leaves, stem, and root tissues, 
carbon loss due to turnover and mortality of vegetation can be 
modeled as described by previous work [43]. 

In order to analyze the outputs simulated by the global coupled 
models, a criteria of regions has been employed by Peng et al. 
(2013; 20 14) [26,2 7] in the regional scale performance (Table 1 
and Figure 1). It has been described by Giorgi et al (2000) [44] as 
follows: (1) the size of the regions vary in the range of a few 
thousand to several thousand km in each direction, so that each 
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Table 1. List of 21 regions. 



ID 


Region 


Abbreviation 


1 


Australia 


AUS 


2 


Amazon basin 


AMZ 


3 


southern Soutli America 


SSA 


4 


Central America 


CAM 


5 


western North America 


WNA 


6 


central North America 


CNA 


7 


eastern North America 


ENA 


8 


Alaska 


ALA 


9 


Greenland 


GRL 


10 


Mediterranean basin 


MED 


11 


northern Europe 


NEU 


12 


western Africa 


WAF 


13 


eastern Africa 


EAR 


14 


southern Africa 


SAP 


15 


Sahara 


SAH 


16 


southeastern Asia 


SEA 


17 


eastern Asia 


EAS 


18 


southern Asia 


SAS 


19 


central Asia 


CAS 


20 


Tibet 


TIB 


21 


northern Asia 


NAS 



doi:l 0.1 371 /journal.pone.0a95282.ta01 

region includes at least a few Earth system models (ESM) grid 
points and thus contains the smallest wavelength of ESM solutions; 
(2) all land areas in the World with a number of regional 
management simple shapes have been approximately covered; (3) 
selection of specific regions is in order to represent different 
climatic regimes and terrain settings. On the basis of such criteria 
of regional selection, we calculate the sensitivities of global and 
regional land carbon storage to climate change and increasing 



atmospheric CO2 concentrations. These sensitivities to direct CO2 
effects and climate change are described as (Pg C per ppm) and 
''iLi (Ps C K ') as the following equations, respectively. The 
calculation method can be referred to a previous study reported by 
Cox et al. (2013) [8]. 



AC 



(1) 



(2) 



Where AC™^ is the change in global or regional land carbon 



storage (in PgC), AT^^^is the change in mean global and regional 
atmospheric temperature (in K) and AC^^'^is the change in 
atmospheric CO2 concentration (in ppm), in response to 
increasing atmospheric CO2 concentrations. AT'£'™°'''is the 
change in mean global and regional atmospheric temperature (in 
K) and \(jChmate change in global and regional land carbon 

storage (in PgC), in response to climate change. These changes in 
temperature, land carbon storage and atmospheric CO2 concen- 
tration averaged over the period of the last 30 years of 130-years 
relative to the period of the first 30 years in all cases of calculation 
for the sensitivities responding the increasing CO2 concentration 
(Yl;) and climate change (Pli). 

In this analysis, R is a correlation coefficient between the annual 
variable (e.g., temperature, precipitation) and the natural sequence 
1, 2, 3,..., n, at a given year [45]. R is calculated by the following 
formula: 



E (xk-x)(k-T) 



(3) 
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Figure 2. Variability in anomalies of annual mean surface air temperature (K) at regional and global scales considering climate 
change alone. 

doi:1 0.1 371 /journal.pone.0095282.g002 



PLOS ONE I www.plosone.org 



4 



April 2014 I Volume 9 | Issue 4 | e95282 



Sensitivity of Carbon Storage to CO2 and Climate 




SSA 



400 
200 
0 

-200 
-400 



' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' 

1S60 ISeO 1900 1920 1940 lOSO 19B0 1860 IBSO 1900 1920 1940 1960 1960 1660 1880 1900 1920 1940 1960 1960 




CAM I 
I ' ' ' I ' ' ' I I I I I 

I860 1660 1900 1920 1940 1960 1980 



WNA 




1860 1860 1900 1020 1940 1960 1980 1660 1880 1900 1920 1940 1960 1960 




... 1 ... 1 ... 1 ... 1 ... 1 ... 1 . 






-i ALA 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




GRL - 



400 
200 
0 

-200 
-400 



1860 1660 1900 1920 1940 1960 1980 
■■■■■■■■■■■■■■■■ 




-HadGEM2-ES 



MPI-ESM-LR 

AVE-MODELS 



I860 1860 1000 1020 1940 1960 1980 

IPSL-CM5A-LR 

CanESM2 



• CESMl-BGC 

BCC-CSMl-1 



Figure 3. Same as Figure 2, but for annual precipitation (mm yr ^). 

doi:10.1371/journal.pone.0095282.g003 
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Figure 4. Same as Figure 2, but for soil moisture (kg m 

doi:10.1371/journal.pone.0095282.g004 
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Figure 5. Variability in atmospheric CO2 concentration (ppm) at rate of 1% per year from pre-industrial values until concentration 
quadruples for a 130-year long simulation, variability in anomalies of global land carbon storage (Cland, gC m"^), soil carbon storage (Csoil, gC 
m"'^), vegetation carbon storage (Cveg, gC m"'^), annual mean temperature (TAS, K), annual precipitation (PR, mm yr^') and soil moisture (iVlRSO, kg 
m"^) only allowing for the direct CO2 effect. 
doi:1 0.1 371 /journal.pone.0095282.g005 



Where xj^ is the annual variable in the time k, n is the sequential 
year, and x is the multi-year mean for this variable; t is equal to 
tlie mean of 1 and n. The linear trend of the variable is presented 



within the period of 1860-1989. When a "significant" linear trend 
for a variable is shown, R must pass the significance level (e.g., 
P<0.05) using the student t-test. 
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Figure 6. Spatial distribution in sensitivity of terrestrial carbon storage to upward atmospheric CO2 concentration simulated from 
six the fifth Coupled Model Intercomparison Project (CMIPS) models, units: gC ppm \ 
doi:1 0.1 371 /journal.pone.0095282.g006 



Results 

The global and regional climate consistently show increasing 
trends in temperature (0.05 °C yr"', R = 0.95, P<0.001), but no 
consensus trends in precipitation (PR) and soil moisture are 
presented at regional scales (Figures 2-4). As shown in Figure 3, 
increases in precipitation in high-latitude regions are shown (e.g., 
in northern Asia with 1.17mm yr~^, R=0.97, P<0.01), while 
decreases in precipitation are located in part of the low latitudes 
such as Amazon basin with -0.50mm yr"^ (R = 0.50, P<0.01) 
and Central America with -1.28 mmyr"^ (R = 0.79, P<0.01). 
These spatial divergences on changes in soil moisture are also 
shown in Figure 4. Only considering climate change, the increases 
in soil moisture is mainly distributed in high latitudes, while the 
area principally in Amazon basin exhibits a considerable decrease 
of soil moisture with -0.23 kg yr"' (R= 0.86, P<0.001). 



Figure 5 shows the interannual variability in modeled global 
annual land carbon storage anomaly, soil carbon storage anomaly, 
vegetation carbon storage anomaly, mean temperature anomaly, 
annual precipitation anomaly and soil moisture anomaly only 
responding to the CO2 increases. As a result, the increase in global 
annual land carbon storage has been presented across the whole 
globe. Such an increase in turn affects the soil carbon storage with 
increasing rate of 20.7 gC yr"' (R = 0.989, P<0.001). 

The geographical distribution of in the direct CO2 impacts on 
the terrestrial and regional land carbon storage is shown in 
Figures 6 and Table 2. Sensitivities of land carbon storage, 
considering the CO2 fertilization eflFect alone, are fairly positive in 
the most terrestrial ecosystem, which means the terrestrial 
biosphere acting as enhanced carbon storage by 1.0 Pg C per 
ppm. The most dramatic increases caused by the increased CO2 
concentrations are mostly located in the regions such as Amazon 
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Table 2. The innpact of increasing atmospheric CO2 concentration on global and regional land carbon storage (units: PgC ppm ^). 



ID 


Acronym 


HadGEM2-ES 


IPSL-CM5A-LR 


CESM1-BGC 


MPI-ESM-LR 


CanESM2 


BCC-CSMI-I 


1 


AUS 


0.04 


0.02 


0.01 


0.1 


0.04 


0.1 


2 


AMZ 


0.22 


0.19 


0.04 


0.24 


0.12 


0.21 


3 


SSA 


0.08 


0.03 


0.01 


0.06 


0.06 


0.06 


4 


CAM 


0.04 


0.02 


0.01 


0.03 


0.02 


0.03 


5 


WNA 


0.04 


0.04 


0 


0.03 


0.02 


0.05 


6 


CNA 


0.04 


0.04 


0.01 


0.02 


0.03 


0.03 


7 


ENA 


0.03 


0.03 


0.01 


0.02 


0.05 


0.03 


8 


ALA 


0.02 


0.01 


0 


0.01 


0 


0.01 


9 


GRL 


0.02 


0.04 


0 


0.02 


0.01 


0.01 


10 


MED 


0.03 


0.03 


0 


0.05 


0.03 


0.03 


11 


NEU 


0.04 


0.05 


0 


0.04 


0.03 


0.03 


12 


WAF 


0.09 


0.08 


0.02 


0.12 


0.12 


0.08 


13 


EAF 


0.09 


0.06 


0.01 


0.14 


0.11 


0.1 


14 


SAF 


0.09 


0.03 


0.01 


0.12 


0.08 


0.06 


15 


SAH 


0 


0 


0 


0.01 


0 


0.01 


16 


SEA 


0.07 


0.13 


0.03 


0.09 


0.13 


0.08 


17 


EAS 


0.09 


0.08 


0.03 


0.06 


0.1 


0.09 


18 


SAS 


0.04 


0.04 


0.02 


0.07 


0.07 


0.05 


19 


CAS 


0.02 


0.01 


0 


0.06 


0.01 


0.04 


20 


TIB 


0.03 


0.01 


0 


0.03 


0 


0.03 


21 


NAS 


0.1 


0.15 


0.01 


0.1 


0.06 


0.09 


22 


WLD 


1.21 


1.08 


0.21 


1.42 


1.09 


1.23 



doi:10.1371/journal.pone.0095282.t002 



basin, northern Asia, western Africa and eastern Africa, which 
have larger vegetation biomass than other regions. However, the 
warming greatly reduces the land carbon storage and the carbon 
sequestration. Excluding Tibet, northern Asia, Alaska and Green- 
land, the impacts of the warming on land carbon storage are 
negative in the terrestrial ecosystems with about —41.6 PgCK~'. 
Additionally, it should be noted that these simulations forced by 
both increasing CO2 concentration and climate change are 
properly diflFerent among six ESMs and regions. The multi-model 
values across the globe range from 0.2 PgC per ppm in CESMl- 
BGC to 1 .4 PgC per ppm in MPI-ESM-LR responding the direct 
CO2 effects and from -17.5 PgC K"' in CESMl-BGC to 
— 58.6 PgC K~' in CanESM2 considering climate change alone 
across the whole terrestrial ecosystem (Tables 2-3 and Figures 6- 
7). ^ 

Figure 7 exhibits the spatial pattern of changes in land carbon 
storage only responding to warming. In response to such a 
warming, the largest loss in land carbon storage is mainly 
distributed in central Amazon basin, southeastern Asia, Central 
America and eastern Africa (Table 3 and Figure 7). Generally, 
negative sensitivities of land carbon storage mosdy located in the 
tropics and Southern Hemisphere. Inversely, areas in northern 
high latitudes primarily exhibit a considerable increase of land 
carbon storage. For example, a significantly enhanced carbon 
storage appears in Tibet with 3.26 gC m~^ yr~' (R = 0.98, 
P<0.001). Increases in land carbon storage in northern Asia and 
Greenland Alaska are also detected (Figure 8). The simulation 
shows a consistent pattern (a negative sensitivity for the land 
carbon storage to temperature across the major land) modeled by 
the major models in terms of the sign of sensitivities of land carbon 



storage to temperature. Spatial discrepancies among regions are 
evidendy presented as mentioned above: increases in land carbon 
storage in the mainly temperature-limited regions, but decrease in 
higher temperature regions. On the other side, in comparison with 
different regions in the Figures 7, it's worth noting that the multi- 
model simulations provide fairly consistent results in most tropical 
areas and Southern Hemisphere (e.g. Australia, Eastern Africa, 
Southern Africa, Western Africa, Southeast Asia, Amazon Basin), 
but inconsistent results in a part of northern latitudes (e.g. some 
show positive values; others show negative values). The inconsis- 
tency is particularly evident in high latitudes of Northern 
Hemisphere, such as Alaska and part of northern Europe. 

Spatially, analysis of responses of land carbon storage to climate 
change also reveals major pattern: negative impacts of rising 
temperature at the low and middle latitudes of the World are 
simulated by the major models in terms of the sign of response of 
land carbon storage; it's positive impacts at the high latitudes is 
presented. The changes are defined as averaged values of the 
period 1960-1989 relative to 1860-1889 in order to calculate 
these sensitivities as mentioned above. In Figure 9, in response to 
rising atmospheric temperature, average zonal-annual (sensitivities 
averaged over the regions and year) patterns and the agreement 
between simulations are deduced: &80% of areas across the whole 
terrestrial ecosystems shows the same sign of land carbon storage 
sensitivities modeled by a majority of simulations (£4/6 simula- 
tions). A majority of the simulations agrees on the negative impact 
of temperature in tropics and positive impact in temperature- 
limited regions in response to climate change. For example, 99.6% 
of Amazon basin and 99.2% of western Africa show the negative 
impacts of increasing temperature on the land carbon storage. In 
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Figure 7. Same as Figure 9, but to the warming (gC m ^ K ^). 
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contrast, the reverse sensitivities to warming are shown botli 
quantitatively and qualitatively in high latitudes of Northern 
Hemisphere. In the most areas, there is litde agreement among 
simulations for the sensitivities to increasing temperature around 
60°N and Greenland. 

Discussion 

In terms of regions (Figures 2-4), the spatial patterns of 
temperature, precipitation and soil moisture are so different. In 
high latitudes and altitudes, temperature and radiation impose a 
complex and varying limitations on vegetation activity [60]. For 
example, in the Western Europe, solar radiation is an evident 
limited factor to the vegetation growth. In the eastern of Tibetan 
plateau temperature change is the main cause to affect the 
vegetation growth and NPP [47]. Accompanied with such 
warming, increases in soil moisture have been reported by Peng 



et al. (2013) [26] in these regions. However, the tropical change 
estimated from the simulations is currendy different from that of 
temperature-limited regions. Indeed water limitation is enhanced 
(e.g. increased droughts) in some of these tropical regions. 
Regionally, the C02-induced climatic change enhances a decrease 
in soil moisture most in Amazon basin and part of arid and semi- 
arid regions (e.g. Mediterranean basin, Sahara and central Asia). 
In some cases, only considering the radiative CO2 forcing, the 
increasing atmospheric CO2 can influence the climate factors (e.g., 
rising temperature). Although the increasing atmospheric CO2 
only influences the climate but not the biogeochemistry, the 
climatic change in turn affects the biogeochemistry. It can be to 
say biogeochemistry can respond to rising temperature and 
changes in other climatic variables. For example, associated 
warming, the increased evapotranspiration produced by increasing 
temperature can also exacerbate drought in Amazon basin [48]. 
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Table 3. The impact of rising temperature on global and regional land carbon storage (units: PgCK"^). 





ID 


Acronym 


HadGEM2-ES 


IPSL-CM5A-LR 


CESM1-BGC 


MPI-ESM-LR 


CanESM2 


BCC-CSMI-I 


1 


AUS 


-1.93 


-0.99 


-0.19 


-6.25 


-3.42 


-3.15 


2 


AMZ 


-13.88 


-9.99 


-7.45 


- 1 2.85 


-10.58 


-8.49 


3 


SSA 


-2.78 


-1.25 


-0.84 


-2.61 


-3.85 


-2.48 


4 


CAM 


-2.04 


-1.53 


-1.49 


-1.9 


-2.01 


-1.52 


5 


WNA 


-1.03 


-0.6 


-0.1 


-0.27 


-1.13 


-0.54 


6 


CNA 


-0.86 


-2.64 


-0.33 


-1.53 


-1.15 


-0.74 


7 


ENA 


-0.57 


-0.79 


-0.29 


-0.49 


-1.62 


-0.83 


8 


ALA 


1.51 


-0.35 


0.01 


0.19 


0.55 


0.61 


9 


GRL 


1.44 


0.39 


-0.02 


0.39 


0.11 


0.33 


10 


MED 


-1.37 


-1.62 


-0.56 


-2.01 


-2.14 


-1.51 


11 


NEU 


-0.59 


-0.35 


-0.19 


-0.76 


-1 


-0.71 


12 


WAF 


-4.38 


-3.14 


-0.97 


-6.17 


-7.61 


-3.82 


13 


EAF 


-4.71 


-0.97 


-0.22 


-6.22 


-4.27 


-3.18 


14 


SAF 


-4.67 


-1.61 


-0.88 


-6.21 


-4.94 


-2.72 


15 


SAH 


-0.03 


0 


-0.01 


-0.37 


-0.07 


-0.24 


16 


SEA 


-5.26 


-5.13 


-1.28 


-3.9 


-6.54 


-3.45 


17 


EAS 


-1.29 


-3.05 


-0.85 


0.45 


-3.74 


-2.93 


18 


SAS 


-0.86 


-1.85 


-1.44 


-3.19 


-3.53 


-0.75 


19 


CAS 


-0.4 


-0.49 


-0.1 


-1.9 


-0.29 


-0.84 


20 


TIB 


1.65 


0.07 


-0.01 


0.67 


-0.3 


0.21 


21 


NAS 


4.94 


0.09 


-0.05 


0.97 


-0.66 


1.24 


22 


WLD 


-36.52 


-35.92 


- 1 7.46 


-53.78 


-58.59 


-35.49 



doi:10.1371/journal.pone.0095282.t003 



As a consequence, enhanced drought can cause reduction in 
captured carbon in hving biomass. 

There is current lack of information on the accurate magnitude 
of the response of terrestrial carbon storage and the affecting 
causes, and thus large diflFerences exist among different ESMs. 
Such differences among models are determined by the used 
approaches (e.g., the models using a biogeochemical approach to 
calculate the terrestrial photosynthesis) [14], which provides an 
indication of the key potential processes controlling the CO2- 
induced and climate-induced carbon uptake /storage [4,22]. For 
example, CESM-BGC simulates the smallest sensitivities of land 
carbon storage to direct CO2 effects. As this model includes the 
CO2 fertilization effect constrained by nitrogen limitation com- 
pared with other models. Other factors potentially influencing 
land carbon storage include warming resulting in the changes in 
land carbon storage, especially at the high latitudes. In the high 
latitudes and altitudes of Northern Hemisphere, the positive 
impact of the increasing temperature on land carbon storage has 
been shown in Figure 7. This is partly due to, in these regions, the 
enhanced growth caused by the elevated temperatures [49]. 
Although the intense warming there increases soil organic matter 
decomposition, soil organic carbon storage has been observed to 
continue to increase [2]. This change is partly caused by the 
increase in vegetation productivity, leading to more turnover 
(litterfall) into the soil [2]. Nevertheless, a qualitative modeled 
intercomparison of changes in total land carbon storage to 
increasing temperature, both increasing vegetation storages and 
soil carbon storage provide useful insights in these temperature- 
limited regions. Reversely, multi-modeled simulations conform- 
ably agree on the negative response of terrestrial ecosystems 



mainly coming from the reduction in carbon storage in tropics and 
most Southern Hemisphere. 

The temperature in this study accounts for a highly important 
role in the land carbon storage for high latitudes of Northern 
Hemisphere. Associated with global warming, changes in carbon 
storage fairly depend on the balance between the input of carbon 
as the net primary production (NPP) and the loss of carbon as 
heterotropic respiration (RH) [46] . Over the past two decades, an 
increase in terrestrial photosynthetic activity has been documented 
across high latitudes of Northern Hemisphere [39,50,51]. NPP has 
an increase with a rate of about 18.4TgC yr ^ (Peng et al 
unpublished data) in association with the warming in the regions 
including Alaska, Greenland, northern Europe, Tibet and 
northern Asia. Thus, accompanied with warming global vegeta- 
tion growth is significantly elevated in these high latitudes of 
Northern Hemisphere. Simultaneously, in these regions covered 
by boreal forest, the lengthened plant growing period has been 
observed (Linderholm, 2006; Zhu et al., 2012). The longer 
growing season somewhat contributes such change, based on both 
the earlier onset of spring [52] and the later ending of autumn 
[50]. Overall, both experimental and modeUing studies suggest 
greening [53] and increased NPP [2,50], which infer a upward 
trend of land carbon storage [2] , responding to global warming. 
This is consistent with our results in high latitudes of Northern 
Hemisphere. In contrast, in the hot regions such as tropics, the 
sensitivities of land carbon storage would be generally negative 
responding to the warming. Such a negative response of carbon 
storage [54] and reduced carbon sequestration [55] have been also 
reported in the previous studies. Indeed, unlike in temperature- 
limited regions, our results suggest that total carbon storage can 
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Figure 9. Six model average of the sensitivities of land carbon storage to climate change and agreement between simulations. 
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not be enhanced by the just rising temperature in current high 
temperature regions (e.g., Amazon basin) (Figures 10-11). This 
may lead to a net carbon source of the terrestrial ecosystems tend 
into under the global warming environment. Such a change may 
be attributed to changes in climatic factors. For example, the 
tropics exhibit higher increase in temperature and decreases in 
precipitation and soil moisture as mentioned above. Consequently, 
reduction in precipitation and soil moisture can result in the 
diminution in carbon storage and thus such negative response is 
presented. This is further evident from the 2005 and 2010 drought 
in Amazon basin [56,57] leading to reduction in land carbon 
storage through decreased vegetation productivity and/or in- 
creased respiration [16]. In addition, heat and drought can 
introduce increase in forest dieback [58] and thus a consequence 
of decreased vegetation production and carbon storage [59]. 
Compared with change in the vegetation storage produced by this 
change in vegetation production, a reduction in soil carbon storage 
can be found due to the decreased litter supply from vegetation 
(leaves, stems and roots) to soil [60]. On the other hand, the 
decreased soil storage is commonly generated by the greatly 
accelerated microbial decomposition with the sequentially rising 
temperature. Particularly, the larger differences in the responses of 
total carbon storage to climate change are presented in these 
regions among the ESMs (e.g., from — 13.9PgCK ' in Had- 
GEM2-ES to -7.5 PgCK"' in CESMl-BGC in Amazon basin), 
which suggests that the estimation of changes in temperature and 
precipitation is especially important in evaluating response of 
carbon cycle of the biosphere. In the same region, the study of 
Cramer et al. (2001) also suggested the importance of climatic 
variables for the captured carbon in biomass [13]. Climatic system 
itself is liable to reduce total carbon storage of terrestrial 
ecosystems. Hence, improvement in knowledge about it is very 
important to reduce and even remove uncertainties in sensitivities 
of land carbon storage at global and regional scales, especially in 
the areas of around 60°N covered by boreal forest. 

Moreover, we compared our results with previous works to 
clarify the magnitudes and spatial pattern of sensitivities of land 
carbon storage to rising atmospheric CO2 concentration and 



climate change at both regional and global scales. The used 
simulations are forced by increasing 1 % yr ' until quadruple CO2 
concentration with/ without the carbon-cycle feedback or radiative 
CO2 forcing. Ranges of multi-model sensitivities commonly faU 
into C MIP model range (sensitivities to increasing CO2 between 
0.2~2.8 PgCppm ' and to rising temperature between 

— 177 20 PgCK~') [19]. The multi-model positive response 

to increasing CO2 concentrations suggests that the carbon 
sequestration is improved owning to CO2 fertihzation effect, 
especially in broad areas of forests. Such an increase can be 
explained by the vegetation increased photosynthesis accompanied 
with the rising atmospheric CO2 concentration. For example, the 
increase of 23.6 gC yr"' (R2 = 0.99, P<0.001) in vegetation 
carbon storage has been shown at the global scale (Figure 5), in 
association with the increased GPP with a rate of about 5.0 g C 
m"2 yr"2(R-' = 0.998, P<0.01) (Peng et al unpublished data). For 
the high latitudes of Northern Hemisphere such as northern 
Europe and at tropical latitudes such as western Africa, which 
contain broad areas of forest, great accumulations of carbon into 
plant biomass appear and the consistent result has be found by 
Cramer et al.(2001), Sitch et al. (2008) and Ito (2005) [13,19,60]. 
In addition, simulated sensitivities of land carbon storage 
responding to climate change have large uncertainties. For 
example, the magnitudes of these sensitivities to temperature in 
CanESM2 and MPI-ESM-LR are significantly larger than that in 
CESM-BGC. Compared with the estimates by Arora et al. (2013) 
[14], the range among multi-model sensitivities to rising CO2 is 
larger, while the range between sensitivities to temperature is 
smaller. Also such an extent to temperature is smaller than the 
range assessed by Friedhngstein et al. (2006) [1] and Zickfeld et al. 
(2011) [23]. Such differences straightforwardly lie in the discrep- 
ancy forcing scenario, the used approach and different but 
plausible representations of the underlying physical and chemical 
processes [1,13,23]. In our study, in all cases sensitivities are 
calculated based on changes from the last 30-yr results relative to 
the first 30-yr results of 1 30-yr simulations. Agreed on simulations 
from both the multi-model differences in dynamic global 
vegetation models (DGVMs)[13,19] and Earth system models 



PLOS ONE I www.plosone.org 



13 



April 2014 I Volume 9 | Issue 4 | e95282 



Sensitivity of Carbon Storage to CO2 and Climate 




-HadGEM2-ES 



ISeO 1880 ISOO 1920 1940 1960 1980 



IPSL-CM5A-LR 

CanESM2 



MPI-ESM-LR 

AVE-MODELS 



Figure 10. Same as Figure 2, but for vegetation carbon storage (Cveg, gC m 

doi:1 0.1 371 /journal.pone.0095282.g01 0 



• CESMl-BGC 

BCC-CSMl-1 



PLOS ONE I www.plosone.org 



14 



April 2014 I Volume 9 | Issue 4 | e95282 



Sensitivity of Carbon Storage to CO2 and Climate 




-HadGEM2-ES 



ISeO 1880 ISOO 1920 1940 1960 1980 



IPSL-CM5A-LR 

CanESM2 



MPI-ESM-LR 

AVE-MODELS 



Figure 11. Same as Figure 2, but for soil carbon storage (Csoil, gC m ^). 

doi:1 0.1 371 /journal.pone.0095282.g01 1 



■ CESMl-BGC 

BCC-CSMl-1 



PLOS ONE I www.plosone.org 



15 



April 2014 I Volume 9 | Issue 4 | e95282 



Sensitivity of Carbon Storage to CO2 and Climate 



(ESMs) [1,14], the changes of land carbon storage across the 
terrestrial ecosystem caused by COg increases of about 563.4 ppm 
account for about 247% of the changes induced by warming of 
about 4.3 K. Generally, the magnitude of responses of NPP versus 
RH to climate change [19,61] is still debated and strength of 
dynamical responses of vegetation growth to increasing CO2 
concentrations is qualitatively different among difiFerent models. 
Such magnitude of NPP has been documented to depend on the 
local water availability [62]. These changes in water avaUabiUty 
depend critically upon uncertain regional aspects of climate 
change projections and are therefore likely to be another dominant 
source of uncertainty [14,19]. Hence the simulated divergences in 
precipitation and soil moisture can partly contribute to the 
variances in land carbon storage from the modeled simulations. 
Overall, these uncertainties can be concluded in the processes: 
temperature dependent RH [6,3] and strength of NPP affected by 
drought and CO2 fertilization [62,64,65], and captured carbon in 
the forest aboveground biomass [66], especially in the Amazon 
forest [67,68]. Hence, it is urgendy needed to pay more attention 
on key processes (e.g. quantifying the strength of COj fertilization 
effect and removing or lessening the uncertainty in climate change) 
and critical regions. For example, for tropical regions, which play 
extremely important role in the terrestrial carbon cycle, great 
uncertainties are even maintained in the competition between the 
direct CO2 effects and climate change due to lacking fully 
understanding driven mechanism in dynamic biosphere compo- 
nent. 

Conclusion 

The results of this study show sensitivities of global and regional 
land carbon storage to rising atmospheric CO2 concentration and 
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